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bji Abstract 

We present the general derivation of the fuh non-perturbative equation that governs the mo- 
mentum evolution of the dynamically generated gluon mass, in the Landau gauge. The entire 
■ construction hinges crucially on the inclusion of longitudinally coupled vertices containing mass- 

JL , less poles of non-perturbative origin, which preserve the form of the fundamental Slavnov- Taylor 

(— 1 1 identities of the theory. The mass equation is obtained from a previously unexplored version of the 

Schwinger-Dyson equation for the gluon propagator, particular to the PT-BFM formalism, which 
involves a reduced number of "two-loop dressed" diagrams, thus simplifying the calculational task 
considerably. The two-loop contributions turn out to be of paramount importance, modifying the 
qualitative features of the full mass equation, and enabling the emergence of physically meaningful 
solutions. Specifically, the resulting homogeneous integral equation is solved numerically, subject 
to certain approximations, for the entire range of physical momenta, yielding positive-definite and 
monotonically decreasing gluon masses. 

PACS numbers: 12.38.Aw, 12.38.Lg, 14.70.Dj 
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I. INTRODUCTION 



The dynamical generation of a momentum-dependent gluon mass jl| has received partic- 
ular attention lately, especially in the light of the important results obtained from a large 
number of lattice simulations. Specifically, the observed infrared (IR) finiteness of the gluon 
propagator and the ghost dressing function (in the Landau gauge), both in SU(2) [j, |3| 
and in SU{3) j4-6|, may be explained in terms of such a nonperturbative mass j7-9|, which 
tames the IR divergences of the Green's functions of the theory (for alternative explanations. 



see, e.g. 



10|). 



The Schwinger-Dyson equations (SDEs) constitute the most natural framework for study- 



ing such a nonperturbative phenomenon in the continuum 



llMlSl]. In particular, the emer- 



gence of IR finite solutions out of the SDE governing the gluon propagator has been studied 
in detail in a series of works (see e.g., 3, H, Q)- Particularly interesting in this context is 
the question of how to isolate the integral equation that determines the momentum evolution 
of the gluon mass [iS^. The main purpose of this work is to present the general derivation 
of the complete gluon mass equation, employing the full SDE of the gluon propagator. 

In order to address this difficult question, we consider the set of modified SDE equations 
obtained within the general framework arising from the synthesis of the the pinch technique 
(PT) |l9-22| with the background field method (BFM) {23], known in the literature 



as the PT-BFM scheme 
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24| . This becomes possible by virtue of set of powerful 



relations, known as the background-quantum identities 25 



These all-order relations 



allow one to express the conventional gluon propagator connecting two quantum (Q) gluons 
in terms of the two additional gluon propagators appearing in the BFM: the propagator 
connecting two background (B) gluons, and the propagator connecting a quantum with a 
background gluon {BB and QB propagators, respectively). This, in turn, permits one to 
use the SDE for the BB or QB propagators, which are expressed in terms of the Feynman 
rules characteristic of the BFM, involving fully dressed vertices that satisfy Abelian-like 
Ward identities (Wis), instead of the typical Slavnov- Taylor identities (STIs). The main 
consequence of this formulation is that the resulting SDE may be suitably truncated, without 
compromising the transversality of the answer. As we will explain in detail in the main body 
of the article, one gains a considerable advantage by considering the version of the identity 
connecting the QQ with the QB propagators, instead of the one connecting the QQ with the 
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BB propagators, used in the literature so far. Specifically, the resulting SDE (to be denoted 
as the "Q-B version" displays the powerful block-wise transversality property known from 
the BB case, but has two additional important features: it contains fewer graphs, and the 
limiting procedure necessary for projecting the result to the Landau gauge is significantly 
less involved. 

Even within this improved framework, one still faces the fundamental question of how 
to disentangle from the SDE of the entire gluon propagator the part that determines the 
evolution of the mass from the part that controls the evolution of the "kinetic" term. In this 
work we present a new unambiguous way for implementing this separation, which exploits 
to the fullest the characteristic structure of a certain type of vertices that are inextricably 
connected with the process of gluon mass generation. 

Specifically, a crucial condition for obtaining out of the SDEs an IR-finite gluon propaga- 
tor, without interfering with the gauge invariance (or the BRST symmetry) of the theory, is 
the existence of a set of special vertices, to be generically denoted by V, that are purely lon- 
gitudinal and contain massless poles, and must be added to the usual (fully-dressed) vertices 
of the theory. The role of these vertices is two-fold. On the one hand, thanks to the massless 
poles they contain, they make possible the emergence of a IR finite solution out of the SDE 
governing the gluon propagator; thus, one invokes essentially a non-Abelian realization of 



the well-known Schwinger mechanism 



27 



28l | . On the other hand, these same poles act like 



composite Nambu-Goldstone excitations, preserving the form of the STIs of the theory in 
the presence of a gluon mass. Recent studies indicate that the strong Yang-Mills dynamics 
can indeed generate longitudinally-coupled composite (bound-state) massless poles, which 
subsequently give rise to the required vertices V [29|. 

It turns out that the very nature of these vertices furnishes a solid guiding principle for 
implementing the aforementioned separation between mass and kinetic terms. In particular, 
their longitudinal structure, coupled to the fact that we work in the Landau gauge, com- 
pletely determines the g^gj, component of the mass equation; this is tantamount to knowing 
the full mass equation, given that the answer is transverse (so, the g^i, part is automatically 
fixed from its g^g^ counterpart). If, instead, one had tried to determine the g^j^^, part first 
(or had taken the trace of the full equation, as was done in [isl]), one would have been 
confronted with a much subtler exercise. Specifically, the cancellation of the quadratic di- 
vergences, implemented by virtue of a characteristic identity, leads to a nontrivial mixing 
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between the g^i, component of the mass and kinetic terms. Even though this separation can 
be eventually carried out, it is significantly more involved and delicate than that of the g^g^ 
components. 

As already mentioned, in the present work we include all fully dressed graphs (one and 
two loops) comprising the corresponding full SDE of the QB propagator. Going beyond 
the "one-loop" dressed analysis is highly nontrivial, because it requires the introduction of a 
new \^-type vertex, never considered before. Specifically, in addition to the ^s related to the 
three-gluon vertices, known from the one-loop case, the V vertex associated with the fully- 
dressed four gluon vertex BQ^ must be included in the corresponding "two-loop dressed" 
diagram. As happens in the one-loop case with the three-gluon V, this new four-gluon 
V vertex must satisfy a very concrete QED-like WI, in order to ensure the transversality 
of the "two-loop dressed" part of the calculation. Interestingly enough, and again as a 
consequence of their longitudinal nature and the Landau gauge, the Wis (and in some case 
the STI) satisfied by all V^-type vertices involved in this problem is all that one needs for 
calculating their effects exactly. This fact clearly constitutes an important simplification 
and bypasses the need to actually construct explicitly the corresponding vertices. 

As a consequence of the novel aspects introduced to our approach, the "one-loop" calcu- 
lation presented in the first part of our derivation f Section IVl) p roceeds in a far more concise 
way compared to the corresponding derivation followed in [l8|, rectifying, in fact, the form 
of the resulting mass equation. The two-loop contribution is considerably more cumbersome 
to obtain, and is expressed in terms of a kernel that, in addition to full gluon propagators, 
involves also the conventional, fully dressed, three gluon vertex {Q^). The two-loop part 
of the mass equation is subsequently simplified by choosing tree-level values for a judicious 
combinations of its ingredients, a fact that allows us to carry out explicitly one of the two 
integrations over virtual momenta. In order to gain insight on the numerical subtleties as- 
sociated with this equation, we first consider its limit at vanishing physical momenta, thus 
converting it into a nonlinear constraint. Already at this level, the contribution from the 
two-loop part appears to be of paramount importance, having far-reaching consequences for 
the behavior of the resulting solutions. The detailed numerical solution of the full equa- 
tion (for arbitrary values of the physical momentum) confirms this impression, revealing the 
existence of positive-definite and monotonically decreasing solutions. 

The article is organized as follows. In Section |ll] we introduce the basic notation, together 
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with some of the most important PT-BFM relations, and explain the derivation and general 
structure of the QB version of the SDE for the gluon propagator. In Section IIIII we present 
a brief reminder of the basic elements appearing in the gauge-invariant generation of a gluon 
mass, with particular emphasis on the role and properties of the V vertices. In Section HVl we 
outline in detail the methodology that allows one to extract from the corresponding SDE two 
separate equations, one for the mass and one for the kinetic term, and explain the advantage 
of selecting out the qfj^q^ component. Section |V] is dedicated to the concise derivation of the 



one-loop version of the mass equation^ 
discrepancy with the result given in 



n addition, we explain in passing the reason for the 
18| : in addition, a brief discussion on the ghost sector 
(only "one-loop" in this new SDE version) is included, explaining why the ghost graphs 
do not affect the mass equation. In Section |Vl] we present the full two-loop calculation, 
organizing the corresponding technical aspects into various self-contained subsections. In 
Section IVIII we present the final form of the integral equation that governs the dynamical 
mass, and discuss some of its general properties. In addition, we calculate an approximate 
expression for the new contribution to the kernel of this integral equation, to be used in 
the ensuing numerical analysis. In Section rVIIII we solve numerically the integral equation, 
determine a family of positive-definite and monotonically decreasing solutions, and study 
their dependence on the value of the strong coupling constant. Finally, our conclusions are 
presented in Section HXl 



II. A SIMPLER VERSION OF THE FUNDAMENTAL SDE 

In a general renormalizable gauge, defined through a linear gauge fixing function of 
the Lorentz type (J-"" = 9^y4^), the all-order gluon propagator AjJ^(g) = 5"'^^^y{q) and its 
inverse read 

A,J(?) = ^ [P,M^-\q') + il/i)qM , (2.1) 
where ^ denotes the gauge fixing parameter (^ = corresponds to the Landau gauge), and 

Ptiv{q) = Qtiv - qi,qvl(i' , (2.2) 

is the dimensionless transverse projector. The scalar form factor A(g^) appearing above 
is related to the all-order gluon self-energy n^y(g); this latter quantity, as a consequence 



P,.(g)A(g2)+e 



qi^qv 
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of the BRST symmetry, is transverse to all orders in perturbation theory, as well as non 
perturbatively, at the level of the corresponding SDE. One has then 

q^U^M) = 0; U^M = U{q^)P,M- (2-3) 

The nonperturbative dynamics of the gluon propagator in the continuum is governed by 
the corresponding SDE. As has been explained in detail in a series of articles, the formula- 
tion of this dynamical equation in the context of the PT-BFM formalism furnishes certain 
distinct advantages over the conventional case; most importantly, it permits the systematic 
truncation of the SDE series, without compromising the BRST symmetry, as captured by 

Within the BFM formalism three types of gluon propagator make their appearance, in 
a natural way: (z) the conventional gluon propagator (two quantum gluons entering, QQ), 
denoted (as above) by A(g^); (m) the background gluon propagator (two background gluons 
entering, BB)^ denoted by A(g^); and {in) the mixed background-quantum gluon propa- 
gator (one background and one quantum gluons entering, BQ), denoted by A(g^). These 
three propagators are related among each other by a set of powerful relations, known as 



Background-Quantum identities 



25 



ism 



30 



26| . obtained within the Batalin-Vilkovisky formal- 



3l|. Specifically, 



A(g2) = [l + G(g2)]2A(g2), 
A(g2) = [l + G(g2)]A(g2), 

A(g2) = [l + G(g2)]A(g2). (2.4) 

The function G(g^), whose role in enforcing these crucial relations is instrumental, is defined 
as the Qf^u form factor of a special two-point function, given by (see Fig. [1]) 

KM = -^^Ca [ A';^{k)D{q - k)H,^{-q,q - k,k) 
Jk 

= 9,uG{q') + ^'L{q'), (2.5) 

where Ca represents the Casimir eigenvalue of the adjoint representation (A^ for SU{N)), 
d = 4 — e is the space-time dimension, and we have introduced the integral measure 



(27r)^ 



d'^k, (2.6) 
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FIG. 1: Definitions and conventions of the auxiliary functions A and H. Tlie color and gauge 



coupling dependence for the field combination shown in the H function is c"'{p)A!'^{r)A1'^{q)^ is gf"''^^. 
White blobs denote connected Green's functions, while gray blobs denote one-particle irreducible 
(with respect to vertical cuts) Schwinger-Dyson kernels. 

with /i the 't Hooft mass. In addition, D°'''{q'^) = 6°'^D{q^) is the ghost propagator, and 
ifj^o- is the gluon-ghost kernel shown. The dressed loop expansion of A and H is shown 
in Fig. [T| notice finally that, in the Landau gauge, an important all-order relation exists, 
which links the form factors G{q^) and L{q^) to the ghost dressing function F{q'^) = q^D{q^), 



namely 



32 



35] 



F-i(g2) = l + G(g2) + L(g2 



(2.7) 



The basic observation put forth in H, 12, 24 1 is that one may use the SDE for A(g^), 



written in terms of the BFM Feynman rules, take advantage of its improved truncation 
properties, and then convert it to an equivalent equation for A(g^) (the propagator simulated 
on the lattice) by means of the first relation in Eq. f l2.4p . The resulting SDE contains a richer 
diagrammatic structure than the conventional one, due to the appearance of a new subset 
of graphs, related to the modified ghost sector characteristic of the BFM; specifically, one 
encounters a new type of ghost vertex, involving two gluons and two ghosts (for the complete 
set of all relevant Feynman rules, see, for example, 12|). 

In this article we present an alternative version of the gluon SDE, which displays all 
known desirable features, and, at the same time, has a reduced diagrammatic complexity. 
Specifically, instead of using A(g^) as the starting point, we consider the SDE satisfied by 
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FIG. 2: The SDE obeyed by the QB gluon propagator. Each box encloses a set of diagrams forming 
a transverse subgroup. Black blobs represents fully dressed 1-PI vertices; the small gray circles 
appearing on the external legs (entering from the right, only!) are used to indicate background 
gluons. 

A(g^), shown in Figj2l and the second (instead of the first) relation in Eq. (12 ■4p . Then, 
the corresponding version of the SDE for the conventional gluon propagator (in the Landau 
gauge) reads 



where the diagrams (a^) are shown in Fig. O 

The crucial points to recognize are (z) one has a reduced set of Feynman diagrams (six 
instead of ten) compared to those appearing in the previous formulation in terms of A(g^) 
(six instead of ten, respectively), and {ii) the most important feature of the PT-BFM 
formalism for our purposes, namely the block-wise transversality imposed at the level of the 
SDE for the gluon self-energy, is still present. In particular, the transversality of the gluon 
self-energy is realized according to the pattern highlighted by the boxes of Fig. [2], namely. 



g^[(ai) + (as)]^. = 0; q'^Has) + (a^)]^, = 0; q^'Ha^) + {ae)U = 0. (2.9) 
III. GLUON MASS GENERATION: A BRIEF REMINDER 



As has been explained in detail in the recent literature, the gauge invariant generation of 
a gluon mass proceeds through the implementation of the well-known Schwinger mechanism 
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in the context of a Yang-Mills theory. This mechanism requires the existence of a very special 
type of nonperturbative vertices, to be generically denoted by V (with appropriate Lorentz 
and color indices] , which, when added to the conventional fully dressed vertices have a triple 
effect: (z) they make possible that the SDE of the gluon propagator yields A^^(O) 7^ 0; 
(a) they guarantee that the Wis and STIs of the theory remain intact, i.e., they maintain 
exactly the same form before and after mass generation; {in) they decouple from on-shell 
amplitudes. These three properties become possible because the special vertices (a) contain 
massless poles and (6) are completely longitudinally coupled, i.e., they satisfy conditions 
such as (for the case of the three-gluon vertex) 

p-'-(q)pf^'f^(r)P'''''{p)V^,^,,,{q,r,p) = 0. (3.1) 



The origin of the aforementioned poles is purely non-perturbative: for sufficient 



y strong 



binding, the mass of certain (colored) bound states may be reduced to zero |36N40| . In 
addition to triggering the Schwinger mechanism, these bound-state poles act as composite, 
longitudinally coupled Nambu-Goldstone bosons, maintaining gauge invariance. This has 
been shown to be realized, within a set of simplified assumption, through a detailed study 



m 



29|. 



For the ensuing analysis it is advantageous to introduce the inverse of the gluon dressing 
function, to be denoted by J{q^). Specifically, in the absence of a gluon mass, we write 

^-\q') = q'j{q% (3.2) 

From the kinematic point of view we will describe the transition from a massless to a massive 
gluon propagator by carrying out the replacement (in Minkowski space) 

A-\q')=q'j{q') /\-^\q') = q' J^{q') - m\q'), (3.3) 

where m'^{q^) is the (momentum dependent) dynamically generated mass, and the subscript 
"m" in Jm indicates that effectively one has now a mass inside the corresponding expressions 
{i.e., in the SDE graphs). 

Gauge invariance requires that the replacement described schematically in Eq. (13.30 be 
accompanied by a simultaneous replacement of all relevant vertices by 

f ^ f' = f^ + V, (3.4) 
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where V must be such that the new vertex T' satisfies the same Wis (or STI) as F, but now 
replacing the gluon propagators appearing on their rhs by massive ones. 

To see how this works with an exphcit example, introducing at the same time some 
necessary ingredients for the analysis that follows, consider the fully dressed vertex BQ^, 
connecting a background gluon with two quantum gluons, to be denoted by Ta^^, With the 
Schwinger mechanism "turned off" , this vertex satisfies the WI 

q''f^^,{q,r,p) = p^J{p^)P^,{p) -r'j{r^)P^,{r), (3.5) 

when contracted with respect to the momentum of the background gluon. The general 
replacement described in (13.41) amounts to introducing the vertex 

Kuui(l^r^P)= "^miq^r^p) + V{q,r,p) ; (3.6) 

then, gauge invariance requires that 

q''V^^u{q,r,p) = m\r^)P^,{r) - m\p')P^,{p), (3.7) 

so that, after turning the Schwinger mechanism on, the corresponding WI satisfied by V 
would read 

q^'KnuiQ^r^P) = (f '^m{q,r,p) + V{q,r,p) 

= [p'^Jmip') - m\p')]P^,{p) - [rV^(r2) - m\r^)]P^,{r) 

= A„i(/)P^,(p) - A~\r')P,,ir) , (3.8) 

which is indeed the identity in Eq. (13.51) . with the aforementioned replacement A^^ — t- A^^ 
enforced. The remaining STIs, triggered when contracting r'^^j^{q,r,p) with respect to the 
other two legs are realized in exactly the same fashion. 

Finally, note that "internal" vertices, namely vertices involving only quantum gluons, 
must be also supplemented by the corresponding V, such that their STIs are also unchanged 
in the presence of masses. To be sure, these vertices do not contain type of poles, but 
rather poles in the virtual momenta; therefore, they cannot contribute directly to the mass 
generating mechanism. However, they must be included anyway, for gauge invariance to 
remain intact. 
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IV. DERIVING THE MASS EQUATION: GENERAL METHODOLOGY 



As explained in the previous section, the special vertices V (and V) enforce the transver- 
sality of the rhs of fl2.8l) in the presence of gluon masses. Specifically, writing the A~^(g^) 
on the Ihs of Eq. (I2.8p in the form given in Eq. (13.31) . one has that 



where the "primes" indicates that (in general) the various fully dressed vertices appear- 
ing inside the corresponding diagrams must be replaced by their primed counterparts, as 
in Eq. (13. 4p . Thus, in graph (ai) the BQ^ vertex will be substituted by F', while in graph 
(ag) both the BQ^ and the type of vertices must contain the corresponding V and V 
components, respectively. In addition, in diagram (05) the primed version of the vertex 
BQ^ will make its appearance. These modifications have an important effect: the blockwise 
transversality property of Eq. (12. 9 p holds also for the "primed" graphs, i.e., when (oj) — )• (aQ. 

The Ihs of Eq. (14. ip involves two unknown quantities, Jm{<f) and m^(g^), which will 
eventually satisfy two separate, but coupled, integral equations, of the generic type 



such that g^/Ci(g^, rn?, A^) — ?■ 0, as — > 0, whereas /C2(g^, m?, A^) 7^ in the same limit, 
precisely because it includes the terms contained inside the V terms. 

In the present work we will focus on the derivation of the closed form of the integral equa- 
tion governing m?{q^). To that end, we must identify all mass related contributions contained 
in the Feynman graphs that comprise the rhs of Eq. (14. ip . Now, with the transversality of 
both sides of Eq. (14. ip guaranteed, it turns out that it is far more economical to derive the 
mass equation by isolating the appropriate cofactors of q^qu/q^ on both sides, instead of the 
g^y, or instead of taking the trace. In particular, to obtain the rhs of the mass equation, one 
must (z) consider the graphs that contain a vertex V , and {ii) isolate the qf^qv/q^ component 
of the contributions coming from the V vertices. 

To explain how one reaches the above conclusion, let us first point out that, clearly, in 
the absence of the V , which contain the massless poles, no gluon mass could be generated; 
so, the gluon mass is inextricably connected with the Vs. In addition, as we will see in detail 



[q'Jm{q^)-m\q^)]P,,{q) 



1 + G(g2) 



(4.1) 




(4.2) 
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in what follows, the longitudinal nature of these latter vertices [viz. Eq. (13. ip ]. coupled to 
the fact that we work in the Landau gauge, force the corresponding contribution to the 
gluon self-energy to be proportional to qi, [see Eqs. fl5.9l) and fl6.13p ]. The only exception 
to this rule is the V that appears inside graph (ag), as part of the vertex F'; however, 
the corresponding contribution is shown to vanish identically in the Landau gauge (see 
Section |VI|) . Thus, if we denote by {aDf^u the V^-related contributions of the corresponding 
diagrams, these latter terms are proportional to q^qu/<f only, namely 

(4).. = (4.3) 

so that 

1 + G(g2) ' ^^■^> 
where the sum includes only the graphs z = 1, 5, 6. 

Similarly, the equation for Jm{q^) may be obtained from the q^q^/q^ component of the 
parts of the graphs that do not contain V components. These graphs are identical to the 
original set (ai) — (ag), but now V — > V„i, A — A^, etc. To avoid notational clutter we 
will use the same letter as before, and the aforementioned changes are understood. These 
contributions may be separated in g^y and q^qu/q^ components, 

{a,),, = g,,AM^) + '^-^B,{q^). (4.5) 

Note that graphs (02) and (04) are proportional to g^y only; so, in the notation introduced 
above, -82(5'^) = B^i^q^) = 0. Then, the corresponding equation for Jm{q^) reads 

qJm[q) i + G(g2) ' ^^'^^ 

with i = 1,3, 5, 6. 

We hasten to emphasize that the fact that we focus on the q^qu/q^ terms, instead of the 
g^j^iy, in no way indicates a potential clash with the transversality of the gluon self-energy, 
which is manifestly preserved throughout. In fact, it is precisely the validity of Eq. (12. Sp 
that allows one to choose freely between the two tensorial structures. The point is that the 
transversality of Eq. (14. ip should not be interpreted to mean that the algebraic origin of 
the terms proportional to g^^u is the same as that of the terms proportional to q^qu/q"^- In 
particular, the qfj,qu/q^ of Jm{q^) and rn?{q^) are easily separable, as the Eqs. (14. 4p and (14. 6 p 
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indicate, whereas their g^i, parts are entangled, and their separation is significantly more 
delicate. 

As a particular example of how the gfj_u part requires an elaborate treatment, while the 
Qfj.Qu/<f does not, let us consider the basic cancellation taking place at the one-loop dressed 
level, enforced by the so-called "seagull identity" : the (quadratic) divergence of the "seagull" 
diagram (02) is annihilated exactly by a very particular contribution coming from graph (ai), 
by virtue of the identity (valid in dimensional regularization) 41 1 



!dA„(/c2 



dfc2 ■ 2 



+ ^ / A„(A;' 



0. 



(4.7) 



Ik '^'^ ^ Jk 

The important fact to realize is that all terms involved in the cancellation enforced by 
Eq. (14. 7p are proportional to g^^, and it is only after they are properly combined that 
their total contribution vanishes (as — )■ 0), by virtue of Eq. (14. 7p : instead, their q^q^jq^ 
counterparts vanish individually, in the same limit. 

In order to appreciate this last point, consider the expression 

(4.8) 



Jk 

where f{k, q) is an arbitrary function that remains finite in the limit q^ — )■ 0. Clearly, 

I^,(q)=g^^A{q') + ^Biq'), 



(4.9) 



and the form factors A{q'^) and B{q^) are given by 



Aiq') 
B{q') 



d — i Jk 
1 



.2 jk ■ qf 



f{k,q), 



d-1 



e-d 



{k-qf 
q2 



f{k,q). 



Then, setting {q ■ k)"^ = q^k'^ cos^ 9, and using that, for any function f{k'^) 

[ cos' ef{k') = ] //(P), 

Jk " Ji^ 

we obtain from Eq. (I4.10p that, as q' — )■ 0, 

1 



^(0) = ^ Jffik); 5(0) = 0. 



(4.10) 



(4.11) 



(4.12) 



Evidently, the function / may be such that the integral defining A{0) diverges, while, for 
the same function, B[0) vanishes. 
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Now, in the context of the "one- loop" dressed calculation, the term I^u{(l) originates from 
graph (ai), with the replacement T — > Tm, A — > A^, as mentioned above. The WI satisfied 
by Tm is that of Eq. (13.51) with J ^ J^, and similar but more complicated STIs hold when 
contracting with respect to the other momenta. Note that what appears in the WI and STIs 
is Jm and not A^- Then, as has been shown in |42], the longitudinal part of may be 
expressed in terms of Jm (and other Green's functions), and as a result, the function f{k) 
(at = 0) is given by 

= -A^(A:^)^^%^- (4.13) 
Now, the point is that, in order to trigger Eq. (14. 7p . /(fc^) should be instead 



,4.4) 



To accomplish this, one adds and subtracts 'w?{k'^) to the f{k) of the ^(0) given in Eq. (I4.12p . 
thus obtaining 

At that point, the first term on the rhs of Eq. (I4.15P goes to the g^^ part of the equation 
of Jmi.'f')] when added to the term (02), which is also assigned (in its entirety) to the g^^ 
part of that same equation, it finally gives rise, by virtue of Eq. (14. 7p . to a contribution 
that is free of quadratic divergences and vanishes in the = limit, as it should. On 
the other hand, the second term on the rhs of Eq. (I4.15P is allotted to the mass equation. 
Thus, unlike B{q^), which unambiguously contributes to the q^q^/q^ part of the equation for 
Jmi.'f') [and satisfies automatically -B(O) = 0], the A{q^) contributes to the g^^ component 
of both equations. 

Let us finally point out that the purely longitudinal nature of the vertices V (and V) 
expressed through conditions such as Eq. (13. ip . combined with the fact that we work in the 
Landau gauge, allows one to obtain all relevant contributions simply from the knowledge of 
the WI (or STI) that these vertices satisfy, without the need to construct them explicitly. 
This is particularly important in the case of the vertex Vypc,r, appearing in graph (05); indeed, 
constructing this vertex explicitly would constitute an arduous task, given the comphcated 
STIs that it satisfies when contracted by the momentum of any of its three quantum legs. 
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(ai) (02) 

FIG. 3: The "one-loop dressed" part of the SDE that contains only gluons. Thick lines represent 
gluon propagators endowed with a momentum-dependent mass. The fully dressed "primed" vertex, 
r', enforces gauge invariance in the presence of such a mass. The symmetry factors are 1/2 and 1 
respectively; we also show for the reader's convenience (in this and the next figures) the color and 
Lorentz indices, as well as the momentum routing used in our calculations. 

V. THE "ONE-LOOP DRESSED" MASS EQUATION: CONCISE DERIVATION 

According to the methodology outlined in the previous section, the one-loop dressed 
contribution to the gluon mass equation stems solely from the V^-part of graph {a'^)^^, to be 
denoted by (a]^)^,^. 

The first simplification stemming from the use of the SDE for the QB propagator (as 



opposed to the BB employed in [18|]) is that the Landau gauge limit ^ = may be taken 
directly, in the part of the calculation related to the masses. Indeed, the only source of 
terms proportional to (which require special care) is the tree-level part of the BQ"^ 
vertex, which only affect the equation for J(g^) (and can be easily dealt with, following the 
procedure explained in p]). Then (see Fig. [3] for the Lorentz and color indices as well as the 
momenta routing used in the following calculation) 

= I 9'Ca I rJ2^(g, -k - q)A"^{k)A^''{k + g)Kp.(g, k, -k - q) , 

where r^*|]^ is the conventional three-gluon vertex of the linear covariant gauges. 



(5.1) 



^^!!lpi(l^k,-k-q) = {q- k),ig^^ + {2k + q)^gap - (2g + fc)„^/3^, (5.2) 



and Aai3{k) is the totally transverse Landau gauge propagator, namely 

Aa^ik) = P^^{k)A{e), (5.3) 

[notice the minus sign difference with respect to our general definition (12. ip ]. Finally, the 
trivial color factor 6°'^ has been factored out. 
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As explained in detail in Section UTTt gauge invariance requires that, when contracted by 
the momentum of the background leg, the vertex Vupa{g, k, —k — q) satisfies the WI of (13 .Tp 
with r = k and p = —{k + q), namely 

g'^Kp.lg, k, -k-q) = m\k)Pp,{k) - m\k + q)Pp,{k + q). (5.4) 

Now, it is relatively straightforward to recognize that {a^)^v is proportional to q^q^/q^ only. 
Indeed, the condition of complete longitudinality of V , given in Eq. (13. ip . becomes 

p..'(^)pap(^)p/3.(^ + g)K'p^(g, fc, -k-q) = 0, (5.5) 

from which follows immediately that 

P'^P{k)P^^{k + g)V;,(g, k,-k-q) = '^ [q^'V^'pAq, k, -k - q)] P^P{k)P^^{k + q). (5.6) 

Thus, interestingly enough, the rhs of Eq. (15. 6p is completely determined from the WI of 
Eq. (13.70 : specifically, using (15. 4p . we get 

P^P{^k)P^''{k + g)i> (g, k,-k-q) = ^ [m\k) - m\k + q)] P''P{k)P^{k + q). (5.7) 

Then, using the elementary tree-level WI 

(I'^tlM^ ^' -k~q) = {k + qfPo.p{k + q)- k'P^^{k), (5.8) 

one can show, after appropriate shifts of the integration variables [i.e., {k + q) ^ k], that 
indeed 

(af ),. = M!i^ /■ ^2(^2) ^ _ ^2] A^''{k)A^p{k + q), (5.9) 

q q J k 

yielding 

= ^ / ^'(^') [ik + q? - k'] A"''(fc)A„,(fc + q). (5.10) 
Q' Jk 

Thus, the one-loop dressed mass equation becomes 



\ [ m\e) [{k + qf - A^P{k)A^p{k + g). (5.11) 
1 Jk 



1 + G(g2) g2 

Notice that at the one-loop dressed level the ghost diagrams (03) and (04) of Fig. |2]should 
also be considered. However, their treatment can be simplified by appealing to some basic 
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properties of the ghost propagator in the Landau gauge, estabhshed through detailed large 



10 



43| . Specifically, we 



volume lattice simulations, as well as a variety of analytic studies 
will take for granted that the ghost propagator D in the Landau gauge remains massless, 
D^^{0) = 0, while its dressing function F is infrared finite, F{0) = c > 0. 

The main implications of these properties for the case at hand is that the corresponding 
fully dressed ghost vertex appearing in graph (03) does not need to be modified by the 
presence of l^-type of vertices. Specifically, in the absence of a gluon mass the vertex Bcc 
appearing in (03) satisfies 

g'^r^ = iD-\k + q) - iD-\k) 

= {k + qfF-\k + q)-eF~\k). (5.12) 

If D{q) remains massless, as we assume, the only effect of the gluon mass is to make the 
dressing function infrared finite, i.e., implement in fl5.12p the replacement F{q) — )■ Fm{q)- 
Thus, for instance, if F{q) ~ Ing^, the gluon mass induces the qualitative change of the 
type Fm{q) ~ ln(g^ + m^), accounting for the aforementioned infrared finiteness of the 
ghost dressing function. The point to realize is that this proceeds without the need to 
modify explicitly, by adding to it a l^-type of vertex; will change only through 
its implicit dependence on the gluon propagators (as well as all other vertices), contained 
inside the diagrams defining its own SDE, which have now become "massive". Given these 
considerations, the result given in Eq. (15. lip exhausts the one-loop dressed case. 

Returning now to Eq. f l5.1ip . we emphasize that this equation differs from the correspond- 
ing equation derived in [18], due to a technical subtlety explained in what follows. Specif- 
ically, the method followed in [l8| consisted in taking the trace of both sides of Eq. (14.11) : 
this in itself is totally legitimate, given the explicit transversality of both sides, but, in the 
light of the comments presented in Section llVt has the disadvantage of mixing the (/^j, and 
the q^q^/q^ components of all contributions. As a result, the parts contributing to Jmi.'f) 
and m?{q^) become interwoven. This, in turn, forces one to resort to an argument analogous 
to that employed in Section HVl namely the judicious "completion" of the seagull identity, 
and the subsequent splitting of Aifi) between Jmi'f') and rn?{q^), but now for nonvanishing 
q^. Specifically, one has that 

Hk,,) ^ - <'= + ^'^;^;f +,^1/^"'*> A(^)A(^ + ,). (5,13) 
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which, of course, in the hmit — )■ reduces to the expression given in Eq. fl4.13p . At 
this point, one must add and subtract an appropriate term that would finally, in the limit 



— )■ 0, trigger again the seagull identity. The term employed in 18| in order to accomplish 
this rearrangement was 

m'^{k + q) -m'^jk) 

instead, the correct term to add and subtract is 

m^(A; + q) — m?{k) 



q2 



(5.15) 



It is relatively straightforward to verify that both choices lead to the same limit as q^ — 0, 
but are obviously different for nonvanishing momentum. In fact, after restoring the in- 
tegration sign and all remaining factors, Eq. f l5.15p furnishes precisely the g^y part of 
Pf^y{q)'m'^{q'^), exactly as expected from transversality, while Eq. (15. 14^ clearly does not. 

After this digression, let us return to Eq. (15. lip . The transition to the Euclidean space 
proceeds by using the standard formulas that allow the conversion of the various Green's 
functions from the Minkowski momentum to the Euclidean g^ = — g^ > 0; specifically 



AE(g^) = -A(-g^); m^(g^) = m^(-g^); GM) = G'(-g^); / = W • (5-16) 

Jk Jk-E 

Then dropping the subscript "E" , we arrive at the final result 

^'(^') = "TTW) ? /^'(^')^"''(^)^"''(^ + ^) [(^ + - • (5-1^) 

Consider finally the limit of the above equation as g^ — )■ 0. Using the general Taylor 
expansion [with y = k"^ and w = {k + g)^] 

f{w) = f{y) + {w- y)f\y) + ^(g^), (5.18) 

(the "prime" denotes differentiation with respect to y, i-e.^fiy) = df{y)/dy), employing 
Eqs. f l2Tp and fli?TT|) . and the fact that, in = 4, L(0) = [3^ , we obtain 



m\0) = -g^CAF{0) I k^ ^\k^)[m\k^)\ . (5.19) 



k 



2 

Note that the result obtained for this limiting case coincides with that found in [18*]. 
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(as) («6) 

FIG. 4: The "two-loop dressed" diagrams. The symmetry factors are 1/6 and 1/2 respectively. 
VI. THE "TWO-LOOP DRESSED" CONTRIBUTIONS 

In this section we will study in detail the "two-loop dressed" diagrams (05) and {a^), 
and their respective contribution to the mass equation (see Fig. H] for the color and Lorentz 
indices as well as for the momentum routing). Note that, in the alternative QB version of 
the SDK equation for the gluon propagator that we consider here, there are no additional 
"two-loop dressed" diagrams. In particular, there are no diagrams involving the fully dressed 
BQcc vertices (the fourths subset in the usual BB version of the SDK), simply because these 
vertices cannot be joined in any way with the conventional tree-level vertices appearing on 
the other side of the (would be) diagram, where the Q-type gluon enters. 

A. General considerations regarding the graph (as) 

Before switching the Schwinger mechanism on, the graph (05), for an arbitrary value of 
the gauge-fixing parameter ^, is given by 

= -'fZr ll^'-'m'^'i^ + k)A^'ik + g)f k + q,-i-k,i), (6.1) 

where the tree-level value of the conventional {Q^) four-gluon vertex is given by 

+ rT"'{g,,g.p-g,o.gM)]- (6-2) 

The fully dressed vertex BQ^ satisfies the following WI (all momenta entering) 

P, t) = ig^ [r'T'T^p^ip, t, g + r) + rV'"''^ p^uit^ r,q + p) 

+ r'"rX.p(r,;^,g + t)], (6.3) 
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where it should be emphasized that on the rhs appear the conventional (and not the BFM) 
trilinear vertices, which satisfy STIs with respect to all their legs, i. e. , 

q"T^,,{q,r,p) = F{q)[p^ J{p^)P:{p)H^,{p,q,r) - J{r^)P^{r)H^,{r,q,p)l (6.4) 



and cyclic permutations 44 1. 



Let us next switch the Schwinger mechanism on. Then, both sides of the WI of Eq. fl6.3p 
must be replaced by "primed" vertices. Now, the "primed" three-gluon vertices appearing 
on the rhs are of the type namely r^^^(g,r,p) = Tmaf,u{q,r,p) + Va^^u{q,r,p), where Tm 
satisfies Eq. (16. 4p with J Jm (and H — )■ Hm, which however we refrain from indicating), 
while V must satisfy, correspondingly, 

g"K^,(g,r,p) = F{q)[m\r')P^{r)H^,{r,q,p) - m\p')P^{p)H^,{p,q,r)], (6.5) 

and cyclic permutations. Note the difference between Eq. (16.51) and the corresponding rela- 
tion satisfied by V, given in Eq. (13. 7p : the latter is an Abelian WI with no reference to the 
ghost sector, while the former is an STI, depending explicitly on the ghost-related quantities 
F and H. 

Then, the only possibility for maintaining the original WI of Eq. (16. 3 p intact is if the 
quadrilinear vertex on its Ihs gets also modified into a vertex satisfying the identity 

?"CS(g, r, p, t) = [r::,X,,(g, r, t) + K^iq, r, p, t) 

= ^g' [r'T'KpM t,q + r) + r^r'X^At. r,q + p) 

+ r'"r'T;,^(r,p,g + t)], (6.6) 

where and V satisfy separately the no-pole (F^) and pole (V^) part of the trilinear 
vertex, respectively. In particular, the V part, which we will be of central importance in 
what follows, satisfies 

Q^'KXiQ, r, p, t) = zg' [r'T^K.^ip, t,q + r) + r^T^X.^t, r,q + p) 

+ r'V"'%.p(r,P,g + t)]. (6.7) 

B. The contribution a^{q'^). 

Let us now focus on the part of the diagram (05) that contains the pole component V 
of the fully-dressed BQ^ vertex F', to be denoted by (a^). The projection to the Landau 
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gauge is straightforward, since there are no ^ ^ terms anywhere in this diagram, and one 
obtains 

i4C = / ^ A"^ W A^'^(^ + fc) A^^(fc + g)V;tT(-g, k + q,-i-kj), (6.8) 

where all gluon propagators assume the transverse form of Eq. (15.31) . 

We can apply the totally longitudinally coupled condition satisfied by V\rap, 

P^{q)P^^{k + q)pP'^{l + k)P^P{l)V)^r.p{-q, k + q,-i-k,i) = 0, (6.9) 

to write (16. 8p . after splitting P^{q), as follows 

(4tt = M°i~% / f^'''m"'i^ + k)A^^ik + q)q'V^-{-q,k + q,-i-k,i). (6.10) 
1 Jk Je 

Using then the WI (16. 3p adapted to the present kinematics, as well as the results 

rbnx rxmr-p{0)amnr _ -"^ ^2^2 jrab/ ^ n n \- 

J J J- /.a/37 -^2^ ^ {9^^'y9af} - 9m9l3l), 

fbmx fxrn-p{0)amnr _ •'^„2^2rafe/ ^ n n \- 

J J J-/,a/37 -^2^ ^ [9^^|39a'y - 9f^-l9al3), (6.iij 

one finds that each one of the terms of the WI gives rise to an integral of the form 

t';{q)=W)q^; J = 1,2,3 (6.12) 
so that Eq. (I6.10p can be written as 

(4e = ^^7'^l5"^^Et.(^^), (6.13) 

with 

W) = \{qp9a,-qa9p-y) f / A"''(£)A'^'^(£ + A;)A^^(A; + g)Kpr(-^ - fc, £, fc); 

q Jk Ji 

t2{q^) = ^{qa9(Sy-q-y9a^) [ [A"P{£)A^'^{e + k)A''^k + q)V,rA^,k + q,-q-e-ky, 
q Jk Je 

W) = \{q,9ap - qp9a,) [ f A"P{i)A^'^{i + k)A'y-{k + q)Vr.p{k + q, -i - k, -q + i). 
q Jk Je 

(6.14) 

Now it turns out that, after the appropriate shifts in the momenta, relabeling the Lorentz 
dummy indices, and applying the Bose symmetry of V^^7, the three terms are actually equal. 
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Then, using the totally longitudinally coupled condition for the vertex V, and the fact that, 
in the Landau gauge, kr/^^'^{k -\- q) — —qr^'^'^{k + q), we get 

= ^%{qa9fi,-qp9a,) f A^-(A; + ?) f A^''{e)A^'^{i + k)^V,,^{-i - kj,k). 
q Jk Je 1^ 

(6.15) 

Then, the integral over £ is a function of k (but not of g), and has two free Lorentz indices, 
a and /3, so that 

/ l^'^p{e)l^P''{l + k)^v^py,{-^-k,^,k)^ A{k'')g'^^ + B{k^)k''k^. (e.ie) 

Jt k 

Therefore, 

t{q') = ^%{qa9fS, - qpQa-y) I A^^ik + q) [A{k')g'^^ + B{e)k"k% (6.17) 

? Jk 

But this term vanishes, regardless of the closed form of A{k'^) and B{k'^), because the 
prefactor is antisymmetric under the exchange a j3 whereas the integral is symmetric. 
Thus, one finally arrives at the important result 

aX{q') = 0, (6.18) 

namely that the graph (05) makes no contribution to the gluon mass equation. 



C. The contributions from graph (ae) 

Let us now consider the graph (og), for a general value of the gauge- fixing parameter ^. 
This graph contains the BQ^ and fully-dressed three-gluon vertices, namely F and V. 
We proceed directly to the massive situation, where these two vertices have been replaced 
by their "primed" counterparts, namely 

X !N^''{^)^f^'^{^ + k)T',^,{-^-k,^,k). (6.19) 

Je 

Evidently, V' and V contain the pole parts V and V, respectively, satisfying the general 
properties mentioned earlier. We will next isolate the terms proportional to V and V , since 
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it is these terms that determine the corresponding contribution of the entire graph (ag) to 
the gluon mass equation. This amounts to writing the product FT' as 

rr = {f^ + v){T^^ + v) 

= f^Tm + VT^ + f^V + VV, (6.20) 
and considering only the last three terms. 

1 . Vanishing of the terms proportional to V 

To proceed with the demonstration, note that if we were in the Landau gauge, i.e., if the 
gluon propagators in Eq. fl6.19p had the fully transverse form of Eq. (15.31) . then V would 
vanish identically, due to its property of complete longitudinality, given that it is an internal 
vertex (Q^-type). The limit ^ = may be taken directly in the part of the graph involving 
the term VV, and therefore this term vanishes immediately. As for the combination T^V, 
one can take directly the limit C, = everywhere, thus making it vanish, except for the term 
that contains the tree-level part of Tm that is proportional to Specifically, the tree-level 
part of the BQ^ vertex is given by 

ri°\(-g, k + q, -k) = r(°\(-g, k + q, -k) - r'Krxi-q, k + q, -k), (6.21) 

where the purely longitudinal "pinch part" is given by 

Kt\{-(1^ k + q, -k) = -gu\{k + q)r - Qrvkx. (6.22) 

The contraction of this term with the propagators A^'^(fc + q)^^-^{k) (with ^ still general) 
yields 

r'Krxi-q. k + q, -k)K^\k + q)A'\k) = ^^AUk + q)+ ^^±^At{k). (6.23) 

In this way, the C,~^ term cancels, making the limit ^ — )■ smooth. So, after setting = 0, 
the second term in Eq. (16. 23^ vanishes, because V will be contracted by three transverse 
projectors; on the other hand, the first term survives, since V is contracted only by two. 
However, as we will see now, this last term finally also vanishes, due to a different rea- 
son. Specifically, denoting this contribution by (a^)/^;^ (thus factoring out the trivial color 
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structure S""^), we have 



\ig'C\{g,^gp,-g,pg^,) J AZ{k + q) lA^^{i)A^^e + k)'^V^,s{~i - k,e,k). 

(6.24) 

Now, the integral contains k but no q, and has two free Lorentz indices, a and /3; therefore, 
it can only be proportional to A{k'^)gai3 and B{k'^)kakp. But, since both these terms are 
symmetric under while the prefactor is antisymmetric, this term vanishes. 

2. The term Oq (q^) 

Let us finally consider the term VTm in Eq. fl6.20p . to be denoted by {0^)^^- It is 
convenient to define the quantity 



Y^'\k)= A''^{e)A^''{e + k)T^,s{-i-kJ,k), (6.25) 
Je 

corresponding to the subdiagram on the upper left corner of (ae). Then, (a^)^iy is given by 

(at),. = \ig'Cl {g,^gp, - g,pg^,) jY^^{k)A'^^{k + q)A'\k)V,rx{-q, k + g, -k) (6.26) 

Then, using once again Eqs. 03.71) . 05.51) and 05. 6p . we obtain 

(agV = -XC\{g,^gp,-g,pg^,)'^^ [ [m'{k) - m\k + q)] Ai{k)A^\k + q)Yf{k) 
^ 1 Jk 



and so 



M^ar(g^) , (6.27) 



4{q') = -ig'C\ {q^gp, - q^g^,) ^ / [m\k) - m\k + q)] A{{k)A^\k+q)Y^^ {k). (6.28) 
4 Q Jk 

At this point it is easy to show that the integral Y is antisymmetric under the a (3 
exchange; thus, given also the antisymmetry of the prefactor under the same exchange, 
one can write 

Y^'{k) = {k^g^ - k^gmik'y, Y^k') = ^1 k^g'.Yfik) , (6.29) 

which gives us the final result 

= [m\k')[ik + qY-e][Yik + q) + Yik)]Aiik)A^ik + q) 

4 Jk 

+ -/-^{q'gs, - 2qsq,) [ m^k^Wik + q) - Y{k)]Aiik)A'"ik + q). (6.30) 
^ q Jk 
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D. Explicit check of the two-loop dressed "blockwise" transversality 



The vanishing of the term may appear somewhat surprising, since, in the PT-BFM 
framework that we employ, it is exactly the V and V type of vertices that allow for the 
appearance of a dynamically generated gluon mass in a gauge-invariant way. Thus, one 
might wonder whether the result (I6.18P is in any way at odds with the characteristic property 
of the "blockwise" transversality, mentioned earlier. 

To show that this is not the case, let us contract diagram (05) with the physical momentum 
q; after carrying out the usual splitting of the BQ^ full vertex T — )■ + V, using the 
result (16.1 8p and applying the WI fl6.3p to the remaining term, we get 



3 
4 



JkJe 

(6.31) 



As far as the contribution of diagram (og) is concerned, we know that it cannot be 
projected directly to the Landau gauge, because the tree-level part of the fully-dressed BQ^ 
vertex contains terms proportional to 1/^. However, proceeding as in subsection I VI C l[ one 
writes 

r_(-g, k + q, -k) = r-s(-g, k + q, -k) - C'Krei-Q. k + q, -k), (6.32) 

where evidently the regular part r^''^^ differs from the usual T' by a tree-level term. For the 
regular term, after using the identity 

q'KlJi-q, k + q, -k) = ^-\k)P,,{k) - A-\k + q)P,,{k + g), (6.33) 

one gets 

9^(«r")J = \^g'cl5'^\g,o.gM-g,pgo.,) [ [[A{k + q) - A{k)]Prs{k)p^'{k + q) 



k Ji 



X A'^^ii)A^^{i + k)T-f{-i-k,i,k). (6.34) 

For the pinch part, after using Eq. (16.230 to cancel the dependence, one obtains (now 
in the Landau gauge) 



2 



Al{k + q)+'-r-^Ai{k) 



X /A"^(^)A^'^(£ + A;)F;^5(-£- A;,£,A;). (6.35) 
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FIG. 5: Diagrammatic realization of the WI for the two-loop dressed diagrams in the Landau 
gauge. One has the cancellations ((15)1 + (a6)i2 = and (a6)ii + (06)21 = 0, while the term (05)22 
vanishes. 



On the other hand, observing that 



+ ^) + + qf 



= A^\k) - A^/^k + q) 

+ [A{k + q)- A{k)]P^''{k)PJ{k + q), (6.36) 



Q l«6 



crpS 
m 



Eq. fl6.35p becomes 

^''^t = ^^/C'l5'^^(^?,/3^?.,-M/37){y^/A](A:)A^(£)A^(£ + fc)r 

- / [A^\k + q)A^^{e)A^-{£ + k)K^s{-e-kJ,k) 

+ / j[A{k + q)-A{k)]P^,{k)P2{k + q)A^^{i)A^^{l + k)T';i:'{ 
Jk Je 



{-i-kj, k) 



k,e,k)j. 

(6.37) 
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Clearly the first term integrates to zero, since, being independent of g, it cannot saturate its 
free index /i, while the third term cancels exactly against f l6.34p . As far as the second term 
is concerned, notice that it still contains a pole part, since the total longitudinality condition 
of Eq. (13. ip cannot be triggered in this case. However, after splitting the full vertex T' one 
finds that the part cancels with the term fl6.3ip . while it is easy to show that the pole 
part V vanishes along the same lines described when dealing with graph . 

The realization of the two-loop dressed "blockwise" transversality in the Landau gauge 
is shown diagrammatically in Fig. [51 



VII. THE FULL MASS EQUATION 

After these rather technical considerations, we are now in position to write down the 
all-order mass equation. Using Eq. fl4.4p . together with the results flS.lOp and f l6.30p . one 
finds 



2 



1 + 
1 + G(g2) g2 



^\e)[[k + qf-e]/\''P{k)/\^,{k + q) 



X {l + lig^CA[Y{k + q) + Y{k)] 



3 g'Cj 1 

4 1 + g2 



{q^gs, - 2qsq,) j m\k^)[Y{k + q) - Af (fc)A^^(A; + q). 

(7.1) 



The transition to the Euclidean momenta can be performed by using the standard formulas 
f l5.16p supplemented with the relation YE^qE) = "iY^—q"^); then one obtains (suppressing 
the "E" subscript as usual) 

m'{q') = - ^ f^^^,^ ^ m\e) [{k + qf - fc^] A"^(fc) A.,(fc + q) 
X |l-^/C^[r(A; + g) + F(fc)]| 

^ ^'^^ ^ {q'gs, - 2qsq,) I m\k^)[Y{k + q) - Y{k)]At{k)A^^{k + q). 

Jk 

(7.2) 



+ 



41 + G(g2)g2' 



Interestingly enough, the full diagrammatic analysis presented in sections |V] and IVIl in 
conjunction with the methodology developed in section [TVl may be pictorially summarized, 
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FIG. 6: Diagrammatic representation of the condensed operations leading to the ah order gluon 
mass equation, where we have introduced the shorthand notation m^(g^) = m'^{q'^)[l + All 
internal propagators are in the Landau gauge. 

in a rather concise way, as shown in Fig O 

Note that the quantities entering in Eq. (17. 2p are bare, and must eventually undergo 
renormalization, following the general methodology outlined in various earlier works [cite]. 
To be sure, the renormalized version of Eq. (17. 2p will involve some of the cutoff-dependent 
renormalization constants Zi introduced during the renormalization procedure, in a way 



analogous to what 



mass (see, e.g. 



45 



lappens in the case of the integral equation for the dynamical quark 



46|): however, for the purposes of the present work, we will ignore these 



constants (replacing them, effectively, by unity, i.e., Zi — )■ 1 ), and will simply consider 
Eq. (17. 2p . assuming that the quantities appearing there are now the renormalized ones. 

As has been explained in section HVt the mass equation derived in Eq. I \7.2h constitutes 
one of the two coupled integral equations that govern simultaneously the dynamics of m?{q^) 
and J{q^) [see, for example, Eq. (14. 2p ]. If the corresponding all order integral equation for 
J(g^) were known, then one could attempt to solve the coupled system, after carrying out the 
additional substitution A(fc^) = [k"^ .J {k"^) +rn? {k"^)]^^ [viz. Eq. (13. 3p ] to all gluon propagators 
appearing inside the various kernels. It turns out that the derivation of the all order integral 
equation for J{q^) is technically far more difficult, mainly due to the presence of the fully- 
dressed four gluon vertex BQ^ [see graph (05) in Figj2], which is a largely unexplored 
quantity, with a complicated Lorentz and color structure, and a vast proliferation of form 



factors. In fact, unlike what happens in the case of the three-gluon vertex BQ^ 



no 



gauge-technique Ansatz exist for this four gluon vertex. Thus, for the rest of this analysis 
[see next section] we will study Eq. (17. 2p in isolation, treating all full propagators appearing 
in it as external quantities, whose form will be determined by resorting to information beyond 



28 



the SDEs, such as the large-volume lattice simulations. Therefore, Eq. (17. 2p is effectively 
converted into a homogeneous linear integral equation for the unknown function m?{q^). 

Evidently, the quantity Y^^{k) introduced in Eq. fl6.25p accounts for the bulk of the 
two-loop contribution, and depends explicitly on the fully dressed three gluon vertex F (of 
the type Q^), in the Landau gauge. This Bose-symmetric vertex satisfies the well-known 
STI (16. 4p and its cyclic permutations, which allow, in turn, for the reconstruction of the 
longitudinal form factors of F in terms of J, F, and the various form factors of the ghost- 
gluon kernel H 4^. Clearly, the inclusion of the (ten) longitudinal vertex form factors into 
Y^^{k), and through it into Eq. (17. 2p . will give rise to rather complicated expressions, whose 
numerical treatment lies beyond the scope of this work. 

In what follows we will simplify this preliminary analysis by considering simply the low- 
est order perturbative expression for F, obtained by substituting tree- level values for all 
quantities appearing in Eq. (I6.25p . and using (16. 29 p . It turns out that even so, the resulting 
expression for Y has sufficient structure to effectively reverse the overall sign of the equation, 
and give rise to physically meaningful solutions for m?{q^). In particular, one has 



Y{e 



1 K 



P<-p^i)pp-^^ + k)Tfpp{-E - fc, ^, k) 



1 k 

1 



\p{(i + ky 



9 



p{(. + ky 



+ kf{i + kY 

(7.3) 



15 / , ^ , 63 
-^7E-log4vr + log-J+- 



(7.4) 



and a straightforward calculation yields (in dimensional regularization, Euclidean space) 

^^^'^ = 3(^ T (; 

where /i is the 't Hooft mass introduced at Eq. (12. 6p . Y{k'^) may be renormalized within the 
MOM scheme, by simply subtracting its value at k"^ = Jp, yielding 

(47r)''4 



(7.5) 



VIII. NUMERICAL ANALYSIS 



In this section we carry out a rather thorough numerical analysis of the mass equation 
derived in the previous sections. 
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To begin with, let us rewrite the equation (for the d = 4 case) in a form that will be 
suited for the ensuing numerical treatment. After setting x = q^, observing that {k + g)^ = 
X + y + 2y/xycos6, and using the measure 

1 



k (2vr)3 J, 



de sm^e I dy y, (8.1) 







we obtain 



^2. ^ ^ _^EM 



/ d^^ sin^^^ / dy y m^{y)^{y)^{x + y + 2^Jxy co'g, 6) 
^ Jo Jo 
X {Aix, y, e)B{x, y, e)[l-C (F (a; + y + 2^^ cos 9) + F (y))] 

- CE{x,y,e) {Y{x + y + 2^ cose) -Y{y))}, (8.2) 



where we have used the approximation 3J] 1 + G{x) ~ F ^{x), and set 



, , ^, X sin^ 9 

A(x,y,9) = 3- 



X + y + 2^Jxyco?, 9 

B{x,y,9) = X + 2y/xy cos9; 

xy + X cos"^ 9{x + 2y/xy cos9) + 2{x + 2y/xy 008 9)"^ 

E{x,y,9) = ■ — ^ . (8.3) 

X + y + 2wxycos 9 



A = (8.4) 



and we have defined (a^ = g /4:Tr) 

OsCa 
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In addition, we have introduced the constant C, multiplying the contribution to the mass 
equation that is of pure two-loop origin. Of course, the value of C corresponding to the 
approximate expression Eq. (17.51) that we employ is fixed, namely 

C = STrCAas] (8.5) 

however, during a significant part of the ensuing analysis we will treat C as a free parameter. 
Thus, essentially one disentangles C from the value of a^, and studies what happens to the 
gluon mass equation when one varies independently and C. The reason for doing this is 
twofold: (z) one has the ability to switch off completely the two-loop corrections (by setting 
C = 0), and (ii) by varying the value of C one may study, in some additional detail, the 
quantitative impact of the two-loop contribution. Specifically, the philosophy underlying 
point {ii) is that, whereas the expression in Eq. (17. 5p furnishes a concrete form for the 
two-loop correction, by no means does it exhaust it; thus, by varying C we basically model, 
in a rather heuristic way, further correction that may be added to the "skeleton" provided 
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by the Y{k'^) of Eq. (17. 5p (for a fixed value of as). Of course, the case where C admits the 
actual value of Eq. (18. 5p will emerge as a special case of this general two-parameter study. 
In view of the ensuing analysis it is convenient to measure C in units of SttCa', to this end, 
we introduce the reduced parameter Cr = C/SttCa, and drop the suffix "r" in what follows. 

The integral equation to solve is a homogenous Fredholm equation of the second kind, 
and can be rewritten schematically as 

m^{x) = -\ d9 dy IC{x,y,9)m\y), (8.6) 
Jo Jo 

where b = oo, but in practice we will choose b ^ 1 but finite. A possible way of solving 
this equation is to expand the unknown function in terms of a suitable function basis, and 



subsequently determine the coefficients of this expansion. In particular |47H49| . using the 
Chebishev polynomials of the first kind T^, one can write 

n 

m\x) = ^ + J2ckTk{x). (8.7) 
k=i 

In order to determine the n + 1 coefficients characterizing the expansion, one discretizes 
X G [0,6] choosing the xj values that corresponds to the extrema of the n^^ Chebishev 
polynomial in that interval, i.e., 

= ^cos (^^(n- j)j + ^, j = 0,...,n. (8.8) 

The problem is then reduced to finding the values of A for which the matrix A + XB is 
singular, where 

A, = ST, ; B,, = 6 j\d j\ /C(x„ y, d) T, , (8.9) 

and 5 = 1, unless z = 0, in which case it is 1/2. Specifically one is looking for the smallest 
positive root As of the generalized characteristic polynomial of the matrices A and —B. 
Provided that A^ exists, one can next determine all the expansion coefficients Ck by simply 
assigning to cq a predetermined value (we choose cq = 1) and then solving the resulting 
reduced system; the corresponding value of the coupling constant can then be obtained 
through Eq. (18. 4p . The solution can be finally rescaled by an arbitrary (positive) constant, 
due to the freedom allowed by the linearity of the equation. 
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A. The one-loop dressed case 



Let us start our analysis from the one-loop dressed case, which corresponds to setting 
C = in Eq. f l8.2p . Specifically, notice that, as x — )■ 0, this equation reduces to the following 
nonlinear constraint 

^'(0) = -^AF(O) dy m\y)lCM; = [^'iv)]' , (8.10) 

where Z{y) = yA{y) is the gluon dressing function. Note that the constraint of Eq. f lS.lOp 
is identical to that derived in 18|; remember, however, that the full one- loop equation for 
arbitrary momenta is different, for the reasons explained in section IVl 

As explained in 18|, the usefulness of Eq. flS.lOp lies in the fact that, already at this level, 
one may recognize the difficulty in obtaining physical solutions {i.e., positive definite in the 
entire momentum range), which can be ultimately traced down to the "wrong" sign in front 
of the equation. In addition, one may explore the conditions that might finally overcome this 
difficulty, without having to solve the integral equation in its full complexity. In particular, it 
is clear from Eq. f lS.lOp that, in order to have a possibility of obtaining physically meaningful 
solutions, the kernel /Ci must display a "sufficiently deep" negative region, which might 
eventually counterbalance the overall minus sign; indeed, we have verified that it would be 
immediate to obtain positive-definite (and monotonically decreasing) solutions if we could 
reverse the overall sign of the equation (or, equivalently, the sign of the kernel). Thus, the 
existence of a negative region in the kernel is a necessary condition for obtaining physical 
solutions; a positive-definite kernel, would exclude immediately such a possibility. Of course, 
as we will see shortly, this condition is far from sufficient. 



Specifically, using the lattice data for A(g^) corresponding to a SU{3) quenched lattice 
simulation ^, one can explicitly verify whether or not (and exactly how) the necessary 
condition described above is satisfied. These lattice data are plotted in Fig. [7] (left panel); 
on the same figure we also plot a fit, whose explicit functional form may be found in various 



recent articles 



18 



35 



461]). On the right panel of Fig. [7] we then show the corresponding 



one- loop dressed kernel /Ci, calculated directly from the lattice data^ and then also using the 



^ The kernel in this case is obtained by first calculating from the raw lattice data {yi,A{yi)) the dressing 
function squared data (yi, Z^{yi)); next the derivative data are calculated as a simple first order central 
difference, i.e., the data plotted are ( ^'+^+^' ^ ^ ^^y^^l-y- ^^''^ )• Finally, errors are calculated through error 
propagation. 
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FIG. 7: {Left panel.) Lattice data for the (quenched) 



H = 4.3 GeV, shown together with the log fit of 



18|, 



35|, 



SU{3) ghion propagator renormahzed at 



461 ] . (Right panel.) The one- loop dressed 



kernel derived from the lattice data, compared to the one obtained from the propagator fit. 

aforementioned fit. One clearly observes the zero crossing of the kernel (at ~ 0.8 GeV^) 
and the corresponding negative region. 

The propagator fit shown in Fig. [7] can be then used to construct the full kernel appearing 
in Eq. (18. 2p . when C = 0; the corresponding full integral equation can be then studied 
by means of the algorithm explained above (for the ghost dressing function F we use a 
continuous interpolator of the corresponding SU{3) lattice data). It turns out, however, 
that the solutions obtained have oscillatory behavior, and display large negative regions, 
ultimately voiding them of any physical meaning. 

Evidently, the negative region furnished by the kernel is not sufficiently deep, or it is not 
located in the optimum momentum region, to counteract the effect of the overall minus sign. 
Clearly, if physical solutions are to be found, the full functional form of the kernel must be 
modified. As we will see in the next subsection, this type of appropriate modification is 
indeed implemented dynamically, when the two-loop corrections are included. 

Finally, in order to avoid any confusion related to the conclusions of this subsection 



and the findings of [18|], let us remind the reader that the equation solved in 18| does not 
coincide with the one solved here; therefore, the (non-monotonic) solutions found in fisl, do 
not correspond to solutions of the present integral equation. 
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FIG. 8: Modification of the sliape of the two-loop dressed kernel K,2{y) with varying C. As the 
latter parameter increases the kernel effectively reverses its sign showing a deep negative well in 
the low momenta region. 

B. The tv^o-loop dressed case: finding physical solutions 

Let us now turn the two-loop dressed contributions back on, by setting C > 0. Consid- 
ering again the x — )■ limit first, one has in this case 



Since, as already mentioned, C will be treated as an independent parameter, one can 
study how the shape of the kernel /C2 changes as C is varied. As can be seen in Fig. |8l 
when C = one is back to the one-loop dressed kernel /Ci of the previous subsection. As 
C increases /C2 displays a less pronounced positive (respectively negative) peak in the small 
(respectively large) momenta region. Next, for C > 0.37 a small negative region starts to 
appear in the IR, which rapidly becomes a deep negative well for y < 0.6, with /C2 becoming 
positive for higher momenta^. Therefore, we see that the addition of the two-loop dressed 
contributions counteracts the effect of the overall minus sign of the integral equation, by 

^ It is also interesting to notice that all the curves meet at a common point , which is determined by the 
general condition 







with ~ 0.67 GeV^. 
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t 0.5 1 1.5 2 2.5 3 

a, = 0.22 

FIG. 9: The curve described by the set of the pairs (C, a^) for which one finds physical solutions 
to the full mass equation (I8.2p . The curve starts from the critical value C ~ 0.56 above which 
exactly one monotonically decreasing solution exists (below C there are no solutions). The value 
C* = as(~ 0.88) corresponds to the case in which Y is kept at its lowest order perturbative value. 
Finally, the gray vertical band represents the value for the quenched strong coupling obtained from 
the 4-loop MOM calculation of 5^] renormalized at /i = 4.3 GeV (a^ = 0.22) with a customary 
10% error; the matching values of C are between 1.55 - 1.65. 

effectively achieving a sign reversal of the kernel (roughly speaking, one has /C2 ~ —JCi)- 
When this analysis is combined with the knowledge gathered from the one-loop dressed case, 
one concludes that there exists a critical value C, such that, if C > C Eq. (18. lip will display 
at least one physical monotonically decreasing solution for a suitable value of the strong 
coupling as- 

In order to see if the picture sketched above is confirmed when x 7^ 0, one can study 
numerically the solutions of the full mass equation (18. 2 p following the algorithm described 
at the beginning of the section (Fig. [9]). Specifically, as shown in Fig. |9l there is a continuous 
curve formed by the pairs (C, a^), for which one finds physical solutions. Indeed, for small 
values of C one has that all eigenvalues A are either negative or complex, and no solution 
exists; this absence of solutions persists until the critical value C ^0.56 is reached, after 
which one finds exactly one monotonically decreasing solution corresponding to the smallest 
positive eigenvalue Xg (with bigger positive eigenvalues giving rise to oscillating non-physical 
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FIG. 10: Typical monotonically decreasing solution of the mass equation (j8.2p . The case shown 
have been obtained for the special values C = 0.88 and 1.6, corresponding value of the couphng 
Os ~ 0.88 and 0.22 respectively. The solutions have been normalized so that at zero they match 
the corresponding (Landau gauge) lattice value A~^(0) f» 0.141 GeV^. 

solutions). However, for values up to C ~ 0.8 the coupling needed to get the corresponding 
running mass is of 0{1), while for the quenched case the expected coupling from the 4-loop 
(momentum subtraction) calculation is as = 0.22 at = 4.3 GeV [50!]. This latter value 
is obtained for C ~ 1.55 - 1.65, whereas for C ~ 0.88 one finds the solution to Eq. fl8.2p 
for the lowest order perturbative value of the coefficient f l8.5p . In general one observes, as 
expected, that as C is increased, as decreases, e.g., for C = 1.1, 1.3 and 1.5 one obtains 
solutions corresponding to the strong coupling values ~ 0.51, 0.35 and 0.26, respectively. 

In Fig. [To] we plot the solutions for the most representative C values, i.e., C = 0.88 and 
C = 1.6 (corresponding to, as already said, as ~ 0.88 and 0.22 respectively); notice that 
we have used the linearity of the equation to normalize the solution in such a way that 
the mass at zero coincides with the IR saturating value found in lattice (Landau gauge) 
quenched simulations 5|], or m^(0) = A^^(O) ~ 0.141 GeV^. As can be readily appreciated, 
the masses obtained display the basic qualitative features expected on general field-theoretic 
considerations and employed in numerous phenomenological studies; in particular, they are 
monotonically decreasing functions of the momentum, and vanish rather rapidly in the 
ultraviolet 111, \l 



51 



52j . It would seem, therefore, that the all-order analysis presented here, 



puts the entire concept of the gluon mass, and a variety of fundamental properties ascribed 
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to it, on a solid first-principle basis. 



IX. CONCLUSIONS 

In the present work we have derived the full non-perturbative integral equation that 
governs the momentum evolution of the dynamical gluon mass. This has been accomplished 
by determining through detailed calculations the mass-related contributions of the "two- 
loop dressed" graphs that appear in the SDE of the gluon propagator, within the PT-BFM 
framework. Thus, the complete equation emerges by including these latter contributions to 
those obtained from the one- loop dressed analysis presented in 18 1. 



We have explained in detail the methodology that allows for a systematic and expe- 
ditious identification of the parts of the SDE that contribute to the mass equation. In 
particular, given the manifest transversality of the full gluon self-energy, we have focused on 
the qf^Qi, component of the answer, thus avoiding a number of technical subtleties related to 
the use of the "seagull identity". Then, the relevant contributions stem entirely from the 
parts of the diagrams that involve the special nonperturbative vertices associated with the 
Schwinger mechanism (generically denoted by V), which are longitudinally coupled, contain 
massless poles, and satisfy relatively simple Wis. In fact, the fully longitudinal nature of 
these vertices, coupled to the use of the Landau gauge, obviates the need of their explicit 
construction; indeed, the desired result may be obtained simply by resorting the the Wis 
that these vertices obey, without the need to invoke their explicit closed form. 

It turns out that the inclusion of two-loop dressed contributions has a profound impact 
on the nature of the mass equation, already at the qualitative level. In fact, the kernel of 
the integral equation acquires an extra term, with respect to the one- loop case [l^, whose 
presence modifies the nature of the obtained solutions. Specifically, even within the simplest 
one-loop approximation, this additional term allows for the appearance of positive-definite 
and monotonically decreasing solutions for the gluon mass m^(g^). A detailed numerical 
investigation of the resulting homogeneous linear integral equation was performed, and the 
dependence of the solutions on the value of the strong coupling as and an appropriately 
introduced parameter C expressing our uncertainty on the two-loop term Ylk"^), has been 
studied. 

Our analysis reveals that there is a critical C value (C ~ 0.56) below which no physically 
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admissible solutions may be found. On the other hand, perfectly acceptable solutions are 
obtained above this value, with the solution corresponding to the strong coupling as = 0.22, 
predicted at /i = 4.3 GeV in the MOM scheme, obtained for C = 1.6. Within the lowest 
order approximation employed for calculating the two-loop dressed term Y, a monotonically 
decreasing mass is obtained for the strong coupling value ~ 0.88. The fact that already 
within this simple approximation for Y one is able to obtain a solution of the full mass 
equation for a relatively modest size strong coupling, should be regarded as a success of the 
entire PT-BFM framework. 

Given the importance of the term F (A;^) for this entire construction, it would be certainly 
important to determine its structure beyond the one-loop approximation used. As already 
mentioned in Section IVIH this can be accomplished, to a certain extent, by making use of 
the Ball and Chiu Ansatz for the three-gluon vertex {mI]. An alternative might be lattice 
itself, given that the three-gluon vertex is (z) computed in the Landau gauge, and (ii) it is 
naturally contracted by three transverse projection operators [see Eqs. ( 16.25^ and (16.28^ . 
Interestingly enough, the resulting structure PPPT is precisely the one that is accessible to 
Landau gauge lattice simulations of the gluon three-point function 53 1. 

The full mass equation, even with the approximate version of y(fc^), provides a natural 
starting point for calculating reliably the effect that the inclusion of light quark flavors might 
have on the form of the gluon propagator in general, and on the gluon mass in particular. 
Specifically, recent studies based on SDEs [M^ as well as lattice simulations for Nj = 2, and 
Nf = 2 + 1 + 1 [551], appear to converge to the conclusion that the quark effects tend to 
suppress the gluon propagator in the intermediate and low momentum region, causing it to 
saturate at a slightly lower point, compared to the saturation point observed in the quenched 
case. Should this picture persists further scrutiny, it would seem to indicate, among other 
things, that the gluon mass increases in the presence of light quarks. It would be therefore 
important to obtain a concrete estimate for the behavior of the gluon mass in the presence 
of dynamical quarks, and especially its value at the origin, i.e., m^(0). We hope to be able 
to report progress in this direction in the near future. 
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